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Abstract 



Fifth order, 0{l3^), exact corrections to the non-singlet electron structure func- 
tion in QED are presented. Calculations were performed in the leading logarithmic 
approximation using the ad hoc exponentiation prescription proposed by Jadach and 
Ward and a recurence formula for the elements of the Jadach- Ward series. A com- 
parison with existing third order, 0{l3^), solutions is also presented. The three next 
elements of the Jadach- Ward series were calculated numerically and parametrized 
with an accuracy better than 5-10~® in the range of x between 0.01 and 1. 
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1 Introduction 



The leading logarithmic non-singlet electron structure function D^^[x,f3) is defined as the 
difference between the densities of virtual electrons and positrons in the initial electron. 
It can be found by solving the Gribov-Lipatov evolution equation [1]: 

= ^(1 - ^) + 1 dvD'^'i^v) ® P{-){x) (1) 



where: 



and 
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P{z) = 6{1 - z) (- + 21n + 0(1 - 6 - z)-±^. (3) 

In equation (2) a[s) is the running coupling constant and s = is the energy scale 
characteristic for the process. In the case of a[s) = a we find that: 

/3(.) = ^ln4. (4) 



The convolution symbol (/) stands for: 

Pl(-) ® P2{-){x) = C C dxidx26{x - XiX2)Pl{xi)P2{x2). (5) 

Jo Jo 

We can rewrite equation (1) in the differential form: 

^-2^1= / dx,dx28{x-x,X2)P{x,)D'''{x2,^) (6) 

op 4 Jo Jo 

with the boundary condition 

D''\x,{)) = 8{l-x) (7) 

which means that for s = the electron does not have any internal structure. More 
information about leading logarithmic calculations of the QED corrections can be found 
in [2]. In reference [2] the perturbative solution to the equations (6)-(7) up to third order 
in (3 has been calculated for various prescriptions of ad hoc exponentiations. 
The main result obtained in this paper is the exact solution to the above equation up to 
the fifth order in (3. The accuracy of the solution is estimated and the comparison with 
the existing third order solutions is presented. I propose also an approximate solution up 
to the eighth order in (3 valid in the range of 0.01< x <1 with accuracy of the order of 
10~^. An interesting symmetry property of the solution is mentioned. 



1 



2 Solution to the GL equation 



According to the prescription of the ad hoc exponentiation procedure proposed by Jadach 
and Ward [3] we look for the solution to the Gribov-Lipatov equation in the following 
form: 



N-l 



where the Gribov function reads: 



n=0 



^ '^^ r(i + /3/2) 2^ ^ 



(8) 



(9) 



and is a solution of equation (6) for a; — )• 1 .At the beginning we eliminate the IR cut-off 
e and arrive at the following evolution equation: 



df3 



+ 



- + -ln(l - a;) - - In J D^^x,fi) 
1 r dy 



I 



A Jx y — X 



1 + ^ D-'{y,^)-2D-'{x,^) 



D^^{x,f3) = DG{x,f3)^x,f3), 



-{l + x')-^x,(3) --lnx^x,(3) 



(10) 
(11) 



Substituting: 
we get for $(a;,/3): 

with the condition 

$(1,/3) = 1. (13) 

We can write the solution of the above equation in terms of a power series in f3. Substi- 
tuting: 

^x,fi) = J2{fi/2rM^) (14) 



y - X 



[3^(i + ^)$(y,/3)-2$(.,/3) 



1 + x' 

1-y 



(12) 



n=0 



one can derive a recurence formula for the coefficient functions of the Jadach- Ward series: 



(j)i{x) 

<t>n+l{x) 



2(1 + -") 

-^[2(1 - a;)' + (1 + 3a;') In a;] 

1 f 1 1 " 1 

< —(1 — a;)A„(a;) — (hnix) In a; H — (1 — a;) —— 

/' in"-*^ (^—^] + '^^^^^ - 2^^ 

Jx y — X \l — xj \ y'^ ) 1 — y 1 — x 
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(15) 
(16) 



(17) 



where 



nl Jx \1 



(18) 



The above recurence formula was first obtained in [4]. A„(a;) can be expressed in terms 
of Nielsen's polilogarithms: 



Ai(a;) = 1 — a; + (1 + x)Li2{l — x) — a; In a; 



and for n> 2 



K{x) = (-ir+i[l - a; + (1 + a;)5„,i(l - a;) + a;5„_i,i(l - a;)] 



(19) 



(20) 



Some of the definitions and formulas concerning the Nielsen's polilogarithms useful in 
following calculations are collected in Appendix A. Using the recurence formula (17) one 
can calculate a few next coefficients of the Jadach-Ward series: 



h{x) 



(1 - xY + -(1 - 4a; + 3x^) In a; + — (1 + 7x^) W a; + (1 - a;^)i:i2(l - x) 



(21) 



1 

96 



<l)i{x) 



384 



5(1 - xY + ^(5 - 24a; + 19x^) In a; + ^(1 - 8a; + 7x^) W x 

+ — (1 + 15a;^) In^ x + (6(1 - x^) + 4(1 + x^) In a;) Li2{l - x) 

+ 12(1 - a;')i:i3(l -x) + 2(1 + 7a;')5i,2(l - a;)] 

7(1 - xY + -(7 - 40a; + 33a;^) In a; + -(3 - 34a; + 31a;^) In^ x 



(22) 



+ — (1 - 16a; + 15a;^) In^ x + —(1 + 31a;^) In^ a; + 8(1 + a;^)i:i^(l - x) 

+ (^14(1 - x^) + 8(1 - xY In a; + ^(3 + 13x^) W x^ Li2{l - x) 

+ (24(1 - x^) + 16(1 + x^) In a;) ^^3(1 - a;) + 48(1 — a;^)iz4(l — x) 

+ (4(1 - 8a; + 7x^) - (1 - 17x^) In a;) ^1,2(1 - x) 

-24(1 - a;2)52,2(l - x) - 5(1 - x')S,,s{l - x)] . (23) 

The coefficient (j)2{x) was first calculated in another way in reference [2]. In order to 
estimate the accuracy of the fifth order solution I have calculated numerically the next 
coefficient of the Jadach-Ward series. The ratio of the fifth to the sixth order perturbative 
solutions of the GL equation using the Jadach-Ward ad hoc exponentiation prescription 
is shown in Fig.l. One can see that the accuracy of the fifth order solution is better than 
1.6 • 10~^ in the hard limit and better than 1 • 10~* in the soft one [x > 0.8). 
In Fig. 2 I compare existing third order solutions, exponentiated by Kuraev-Fadin (KF) 
[8] (for the explicit third order result see [2]) and Jadach-Ward (JW) prescriptions. Both 
solutions are normalized to the fifth order solution according to the JW exponentiation 
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procedure. The comparison shows that the JW prescription is closer to the exact result 
than the KF one, especially in the soft limit. 
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Figure 1: The ratio of the fifth (-D^^g) and sixth (-D^^g) order exponentiated solutions 
of the Jadach-Ward (JW) type for the LL non singlet electron structure function. The 
calculation was done for f3 ~ 0.11 i.e. ^/s = 92GeV for the LEP case. 
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Figure 2: The ratio of the third order approximate solutions of the KF (open circles) 
and the JW (full circles) types to the fifth order solution of the JW type. The calculation 
was done for /3 ~ 0.11 i.e. ^/s = 92GeV for the LEP case. 



In Fig. 3 I present all the calculated coefficients of the Jadach-Ward series and the 
next three (^5, 4'7) which I have determined numerically. I have parametrized these 
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last three coefficients using the following function: 

\ pi(6) + pi(7)x + pi[8)x^ + pilyjx'^ + Pi[lu)x^ J 

(24) 

Parameters pi are collected in Table 1. The accuracy of the functions (j)i{x) for i =5, 6, 
7, is better than 5-10~^, in the range 0.01 < a; < 1. 



Table 1: Numerical values of the parameters pi in equation (24). 



i 


5 


6 


7 




(1) 


0.641 


664 


028 


0.003 


843 


748 


-0.000 


103 


279 


Pt 


(2) 


-0.619 


988 


501 


0.006 


263 


836 


-0.002 


208 


024 


Pt 


(3) 


0.224 


513 


784 


-0.033 


976 


652 


-0.061 


516 


549 


Pt 


(4) 


-1.589 


189 


529 


-0.210 


926 


712 


0.130 


086 


362 


Pt 


(5) 


-2.087 


939 


501 


0.156 


033 


173 


0.759 


156 


644 


Pt 


(6) 


0.227 


220 


744 


0.010 


454 


907 


0.000 


757 


690 


Pt 


(7) 


110.025 


192 


261 


-2.389 


885 


902 


-0.249 


451 


116 


Pt 


(8) 


68.255 


920 


410 


-35.333 


641 


052 


-7.883 


694 


649 


Pt 


(9) 


260.508 


605 


957 


-30.422 


819 


138 


-264.039 


489 


746 


Pt 


(10) 


195.836 


334 


229 


-204.567 


718 


506 


-1080.647 


705 


078 


Pt 


(11) 


-0.005 


347 


746 


0.000 


187 


368 


-0.000 


658 


672 



The shape of the approximate function reflects the fact that the fraction of the subse- 
quent coefficients is very close to —1 (the following coefficients are of the opposite sign) 
in a wide a;-range and the singularity of the z-th function for a; — )• is of one power of 
the logarithm stronger than the (z— l)-th one. 

From Fig. 3 we see that the coefficients of the JW series appear to converge to two dif- 
ferent functions: one for the odd and another one for the even coefficients. These limits, 
however, have an interesting property: they are very symmetric with respect to the x axis 
{(j)i{x) = 0). The a;-range [xq < a; < 1) where the symmetry follows with good accuracy 
increases with increasing coefficient function index. 

In order to test how close the subsequent coefficient functions of the JW series are we 
present in Fig. 4 the modules of their fractions. One can see that the coefficients ^5,^6 
and ^7 are very close to each other in the soft a;-range (for x > 0.5 the difference is of the 
order of 0.01%, and for x > 0.01 the difference is up to a few percent). 
The next coefficient, is not drawn in Fig. 3 and Fig. 4 because the difference between 
them and ^7 is smaller than 0.1% in the hard range and smaller than 0.01% in the soft 
one and would be not seen in the figures. 
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Figure 3: Comparison of all obtained (analytically and numerically) coefficients of the 
JW series. 
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Figure 4: Comparison of the fractions of modules of the coefficients of the JW series. 
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3 Conclusions 



Summarizing, we now have analytical, fifth order, corrections to the non-singlet 

electron LL structure function. The result was obtained according to the ad hoc expo- 
nentiation prescription of the Jadach-Ward type. The accuracy of the new solution is of 
the order of 10~^ in the hard limit and of 10~* in the soft limit. The comparison with the 
existing third order solutions shows that the Jadach-Ward exponentiation provides an 
especially good aproximation. An approximate solution up to the eight order in f3 (with 
numerical accuracy of the order of 10~^ in the range 0.01 < a; < 1) was also proposed. 
The higher order coefficient functions seem to be very close to ^7 in modulus especially 
in the soft limit. The subsequent coefficients are of the opposite sign in the wide x-range 
and seem to converge to two very symmetricUy placed limit functions: one for the odd 
and another one for the even coefficients. 
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Appendix A 



Generalized Nielsen's polilogarithms are defined as 



Sn,m ( 2^ ) 



(n-1) 



!m! 



Jo t 



(A.l) 



Sn,m{y)i defined only for positive integers n and m, is real for real y < 1. From the 
definition of Sn,m{y) one can find its derivative and integral: 



-^Sn,m{y) = -Sn-l,m{y) 

ay y 

/ 'S'n,m(2^) — 'S'n+l,m(y) 

Jo X 



In particular: 



-^Li2{y) 
dy 



d 
dy 



SiAy) 



ln(l - y) 
y 

Liijy) 
y 

ln^(i - y) 
2y 



(A.2 
(A.3 

(A.4 
(A.5 
(A.6 



where Lin{y) = S'„_i,i(y). Some of relations between polylogarithms of different argu 
ments: 



Li2{l -y) = -Li2{y) - ln(y) ln(l - y) + C(2) 

Li2{--^) = -Li2{y)-lln\l-y) 
l-y 2 

Lisil -y) = -Si,2{y) - ln(l - y)Li2{y) 



--ln(y)V(l-y) + C(2)ln(l-y) + C(3) 



(A.7 
(A.8 

(A.9 



y 



Lis{--^) = S,,2{y) - Lis{y) + ln(l - y)Li2{y) + -ln^(l - y) (A.IO 
l-y b 



SiM -y) = -Lisiy) + Hy)Li2{y) + - ln(l - y) \n\y) + C(3) 



^i,2(- 

where C(^) = Lin{l) 



j^) = S,,2{y)--ln'{l-y) 
l-y b 



(A. 11 
(A. 12 
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Some of definite integrals used in the calculations: 

/'Un(c + ey)^ _1 ^ /fee - ae V / a + b\ ^/ « + ^\^^-/ «e \ 
^0 a -\- by b \ b J \ a J \ ae — be J \ae — be J 

(A.13) 

r Mi-y)Hi--y) ,^^SM + Li,{a) (a.14) 

Jo y 

C MylMinM^y = 25i,2(a) - Li^{a) + ln(l - a)[Li,{a) - C(2)] (A.15) 
^0 l-y 

I dy = -5i,2(l - X)- Li3{l - X) (A. 16) 

Jx y — X 

f IMLiytilMl^,y = - .) - - (A.17) 

Jx y ~ X 

S.,2{-^ - y) - S,,,{1 - x) ^ _ ^ ^^^^^^ _ _ 1^.2^^ _ ^^_^g^ 
7a; y — a; z 



More detailed information about polylogarithms and associated integrals can be found in 
[5],[6]. 
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